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Abstract 

We simulate the fragmentation processes in the Ca + Ca collisions at the bom¬ 
barding energy 1.05 GeV/u using the Lorentz covariant RQMD and the non-covariant 
QMD approaches, incorporated with the statistical decay model. By comparing the 
results of RQMD with those of QMD, we examine the relativistic effects and find that 
the multiplicity of the a particle after the statistical decay process is sensitive to the 
relativistic effects. It is shown that the Lorentz covariant approach is necessary to 
analyze the fragmentation process even at the energy around Lii a b= 1 GeV/u as long 
as we are concerned with the final observables of the mass distribution, particularly, 
the light fragments around A = 3 ~ 4. 
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The main aim of the high and intermediate energy heavy-ion physics is to investigate 
properties of nuclear matter under extreme conditions, especially to determine the nuclear 
equation of state (EOS). The numerical simulation approaches, such as BUU |], |2| and 
QMD ||, are most promising in describing the time-evolution of the complex system. 
There are two important ingredients of transport theories: a mean-field for nucleons and 
in-medium baryon-baryon cross-sections that accounts for elastic and inelastic channels. 
By searching the appropriate parameters of the mean-field which can fit the experimental 
data, one expects to know the nuclear matter properties 

Within the framework of BUU-simulations we have succeeded to predict/reproduce 
particle production data in heavy-ion collisions and to clarify their reaction processes |5]. 
In spite of these successes it is not always possible to extract the nuclear EOS from the 
results of theoretical calculations without ambiguities for other model-inputs. 


In order to achieve the highly dense nuclear matter, one needs to make experiments 
of nucleus-nucleus collisions above about 1 GeV/u. In such a relativistic energy region, 
there are important relativistic effects other than kinematics. The first author || has 
first suggested that the transverse flow is influenced by the Lorentz contraction of the 
initial phase-space distribution and the Lorentz covariance of the interaction. In further 
developments the Tuebingen group have made more detailed discussions on relativistic 
effects in the dynamical evolution of the nucleus-nucleus collisions, comparing the results 
of the covariant and non-covariant QMD [j?. Q. 


As for the BUU approach, the full relativistic formalism, so called the Relativistic 
BUU (RBUU) 0, has been constructed by incorporating the relativistic a-co model m 
with the BUU. The first author have introduced the momentum-dependent mean-fields 
into the Relativistic BUU approach and succeeded to reproduce the experimental data 
of the transverse flow and the subthreshold K + - production around £’i a b= 1 GeV/u 
simultaneously |TT|. However the fragmentation has not been well discussed at this high 
energy, particularly from the view point of the Lorentz covariant effects. 


Recently, as a new round of the high energy heavy-ion physics, fragment distributions 
are investigated most vigorously by a lot of experimental groups |T2| . The QMD approach 
is the most useful one for this study. In the low energy regime around several tens MeV/u, 
we can nicely reproduce experimental data by simulating the statistical decay [|13[] from 
excited fragments obtained at the end of the QMD calculation |T3. For the high energy, 


however, we have the following problems when we introduce the relativistic kinematics 
and the Lorentz transform into the non-covariant QMD approach. First, the increase of 
the initial density due to the Lorentz contraction make an additional repulsion through 
the density-dependent force, which causes the spurious excitation and the unphysical 
instability of initial nuclei. It is a serious problem at the ultrarelativistic energy 0 . 
Second we cannot correctly evaluate the internal energies of fast-moving fragments at 
the end of the QMD calculation, because the non-relativistic mean-field used in QMD 
is variant under the Lorentz transformation. Hence a fully Lorentz covariant transport 
approach is desired in the relativistic energy region. 


The Relativistic QMD (RQMD) approach 0 [T7j|| must be the most useful theoretical 
model for this purpose; it is formulated to describe the interacting IV-body system in 
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a fully Lorentz covariant way based on the Poincare-Invariant Constrained Hamiltonian 

0 


Dynamics 


The position and momentum coordinate, g* p,-, of the i-th nucleon are 


dehned as four-dimensional dynamical variables and the functions of the time evolution 
parameter r. The on-mass-shell constraints are given by 


H, = pj - rri\ - Vi = 0, 


( 1 ) 


where m, and Vi(qj,Pj) are a mass and a Lorentz scalar quasi potential. The detailed form 
of the quasi-potential is determined by the requirement that it is corresponding to the 
non-relativistic mean-field in the low energy limit |]T5], [16|. Whereas in the non-relativistic 


framework the argument of the potential is a square of the relative distance between two 
nucleons r?, in RQMD we take it as a square of the relative distance in the rest frame of 
their CM system as 

(qh ■ Pij f 


Qth ~ ~Qij + 


Pij 


<li.i <P - <lj- Pij=Pi + Pj- 


( 2 ) 


(3) 


The change from r? 


to 


-?Ti 


L . causes the direction dependent forces and the Lorentz 
contracted phase-space distribution of fast moving matter. In this formulation the time 
coordinate q { - is distinguished from the time evolution parameter r. From that one can 
take into account the retardation effect approximately for the short range interaction such 
as usual nuclear force. Using the above on-mass-shell and the time-fixation constraints 
00- we can describe the Lorentz covariant motions of nucleons. These constraints are 
chosen to be completely consistent in the non-relativistic framework at the non-relativistic 
limit (rrii —► oo). 


In this paper, we calculate the fragment distribution in the Ca + Ca collisions at 
the bombarding energy E\ ah = 1.05 GeV/u with the RQMD approach plus the statistical 
decay model, and then discuss the Lorentz covariant effects in this reaction by comparing 
the results with those of the non-covariant QMD approach. In order to compare these 
two models precisely, we have prepared an integrated code which includes the RQMD and 
QMD parts. The common parts of this code, for example, the prescription of the wave 
packets, the initialization of the ground state, the two-body effective interaction and the 
collision prescription are shared by the two models. By this code we can minimize the 
differences between the RQMD and QMD calculations caused by any accidental difference 
in the coding. Thus the differences observed in the final results are regarded as an evidence 
of the Lorentz covariant effects or the difference of the boosting method discussed below. 
The details of the code and the way of the actual calculations are as in the following. 


As for the basic part of the RQMD code, we follow the prescription in Ref. 0 


while for the QMD and the statistical decay codes the basic prescriptions are taken from 
Ref. [0, 0. For the present purpose, however, we use the relativistic kinematics even in 
the non-covariant QMD approach. Starting from the initial distributions made with the 
cooling method [181, we boost them according to the bombarding energy. Then we perform 
the RQMD and QMD calculations and obtain the dynamical fragment distribution. We 
boost each dynamical fragment to its rest frame and evaluate its excitation energy in the 
way of Ref. [0, which is used in the statistical decay calculation. We used the predictor- 
corrector method to integrate the equation of motion for both models to retain strictly the 


3 
























energy conservation, which is important to estimate the excitation energy of the fragments. 
For the two-body effective interaction, we use a Skyrme-type interaction with HARD EOS 
(the incompressibility K = 380 MeV) and the symmetry force (the symmetry energy is 25 
MeV) [0, while we omit the Coulomb force for simplicity. The widths of the wave packets 

lH is 


are taken from the values for Ca + Ca in Ref. |2C|. The Cugnon parametrization 
used as a cross-section of two baryon collisions for elastic and inelastic channels with the 
frozen-delta approximation. The delta-mass distribution and the angular distribution of 
the inelastic channels are given in Ref. . 


As explained before, we need to transfer the phase-space distribution twice between 
the rest and moving frames. The Erst boost is that of the initial nuclei, namely the initial 
boost, to give a bombarding energy, and the second is that of the moving fragments 
at the end of the dynamical calculation into their rest frames, namely the final boost, to 
evaluate their excitation energies for the statistical decay calculation. In the non-covariant 
QMD a way of these boosts is completely ambiguous; we cannot conclude whether the 
boost should be of the Lorentz or of the Galilei. (Note that the Galilei boost does not 
change relative momenta and positions of nucleons inside the nuclei and fragments in any 
frame. The total momenta of the nuclei and fragments are determined by the Lorentz 
transformation even in this case.) In fact we do not have any definite way to separate 
a total energy of a fragment into its translation energy and internal energy in the non¬ 
covariant formalism using the relativistic kinematics and the non-relativistic mean-held. 
In addition the non-relativistic interaction prefers the spherical density distribution, and 
then even a fast-moving nucleus or a fragment cannot hold a Lorentz contracted shape. 


To examine the relativistic effects in the simulation we compare two kinds of boosts for 
the non-covariant QMD calculation as follows. One (QMD/G) is defined by the method 
that both the initial and the final boosts are of Galilei. The other way (QMD/L) is that 
the both boosts are of Lorentz. The actual equations are given as follows. For the initial 
boost 


r; = a ■ r? + R, 

Pi = b- p° + P, (4) 


where r? and p? are the initial position and momentum of the i-tli nucleon in the rest 
frame, and R and P are the mean position and momentum of the initial nucleus, where 
the initial mean momentum is evaluated from the bombarding energy with the Lorentz 
transformation. In the Galilei boost the factors a and b equal to the unit matrix, while in 
the Lorentz boost they are given as 


a = b 1 


7 


1 + (7 — 
1 

\/l — v 2 


(l 0 0 \ 

0 10 , 
VO 0 7/ 


(5) 


with the boosting velocity v and the unit vector e z along the beam direction. For the final 
boost 


frag 

r= 



4 














p['“ = b, ■ (Pi - P'“ s ), 


( 6 ) 


where R frag and P frag are given as 

Rfrag = £ 

i G frag 

P"‘ e = E PiAW 


(7) 


iGfrag 


Here the factors Of and 6f equal to the unit matrix in the Galilei boost, whereas they are 
given in the Lorentz boost as 

Of = ^ f _1 = 1 + (7f - l)e„ f t e„ f , 

1 

7f = ; 

1 - v\ 


( 8 ) 


with the fragment velocity Vf and the unit vector e 1)f along the fragment velocity. 


In RQMD the definition of the fragments produced in the dynamical process is not 
trivial due to the multi-time description. We thus approximately define the fragments 
at the end of the RQMD simulation in the following way. First, we search candidates 
of the dynamical fragments, where all nucleons have neighbours with the condition that 
—qfj < 16 fm 2 . We then boost the nucleons inside each candidate into its rest frame by the 
Lorentz transformation. With the above condition, in addition to the cluster separability 
condition employed in RQMD |15|, we can roughly define the cluster in which the time 


competent of q\ of all nucleons are not so much different each other. Finally we again 
search dynamical fragments in this frame with the same procedure used in the non¬ 
covariant QMD |2(J. We note that the final result of the mass distribution of the fragments 
obtained with this method does not change within the statistical errors, if we vary the 
value of the above condition from 14 fm 2 up to 18 fm 2 . We also apply this method to the 
non-covariant QMD with the two kinds of the boosts mentioned above. 


Before studying relativistic effects, we check the low energy limit where RQMD and 
QMD should give the same result. We show in Fig. 1 the impact-parameter dependence 
of the a multiplicity in Ca + Ca collisions at the bombarding energy fq ab = 50 MeV/u. 
This figure surely shows that the results of RQMD agree with those of QMD at the low 
energy. 

Next we simulate Ca + Ca collisions at the bombarding energy E\^— 1.05 GeV/u. In 
the left-hand side of Fig. 2 we show the impact-parameter dependence of the dynamical 
single proton (a) and the a (b) multiplicities obtained in RQMD, QMD/G and QMD/L 
calculations before the statistical decay. It can be easily seen that there is no clear differ¬ 
ence between the results of these calculations. Though the Lorentz contraction in QMD/L 
causes the unphysical instability in the initial distribution as mentioned before, the figure 
shows that this instability is not so significant as to emit nucleons in the dynamical stage 
at this energy around -Fi a b= 1 GeV/u. 

In the right-hand side of Fig. 2, on the other hand, we show the same quantity as 
in the left-hand side but after the statistical decay. One can see obvious differences in 
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the results of these models. The multiplicity of single protons (Fig. 2c) decreases with 
increasing the impact parameter b and becomes almost zero at b « 8 fm in RQMD and 
QMD/G, while it has still a finite value in QMD/L. In such peripheral region two-baryon 
collisions rarely occur so that this result of the QMD/L simulations is unphysical. This 
unphysical behavior is much enhanced in the a multiplicity (Fig. 2d). These are due to 
the spurious excitation in the QMD/L calculation caused by the initial Lorentz boost and 
the non-covariant treatment of the interaction. The QMD/G approach seems to be free, 
at least, from this problem. 


In Fig. 2d the multiplicity of the a particle is shown as its impact-parameter depen¬ 
dence. By comparing the results between before (b) and after (d) the statistical decay 
process, it is clear that the a particles are produced dominantly in the evaporation pro¬ 
cess from the excited fragments. These excited fragments are coming from the spectator 
zone, since the participant breaks mainly into single nucleons in this energy region. Thus 
the a multiplicity is large in the middle impact parameter region, as seen in Fig. 2d. In 
this middle impact parameter region, around b = 4 fm, we can see the clear difference 
between the results of RQMD and QMD/G, which amounts to roughly 50 %. 


In RQMD, the Lorentz contracted initial distribution compacts the size of the partic- 

[16 


and nucleon momenta in this zone, 
and makes a stronger repulsion 


ipant zone and increases the maximum density 
Then the participant in RQMD is more heated [^3 
than in QMD/G. Furthermore, the reaction time is shorter in RQMD |Rj||, and the spec¬ 
tators and participant are separated more clearly in the momentum space in RQMD than 
in QMD ||. As a result, the average temperature of the spectators is lower in RQMD 
than in QMD/G. This difference between RQMD and QMD/G also affects the proton 
multiplicity in Fig. 2c. However, the proton is mainly produced in the dynamical process 
(cf. Fig. 2 a and c). Thus this relativistic effects appear only in the a multiplicity. 


Finally, we show the total cross-section of fragments with the mass number A after 
the statistical decay in Fig. 3. Around A = 3 and 4, which are mainly evaporated from the 
excited spectators, we can see the clear difference about factor of two between the results 
of RQMD and QMD/G. It is noted that before the statistical decay the mass distributions 
of the fragments are almost the same in RQMD and QMD/G within the statistical errors. 
Therefore, the difference around A = 3 and 4 is attributed to the difference of the internal 
energies of the fragments produced by these models. 


In summary, we calculate the fragment distribution using the RQMD and QMD ap¬ 
proaches combined with the statistical decay model. The mass distributions of the frag¬ 
ments are almost equivalent in the dynamical stage between RQMD and QMD at the 
energy around 1 GeV/u. This agrees with the conclusion of Ref. [||], where it is claimed 
that the relativistic effects on the observables like the resonance production, density and 
rapidity distributions are less significant at the energy f? lab < 2 GeV/u. However, the dif¬ 
ference appears in the internal energies of the fragments produced by these models. We 
have shown that the multiplicity of the a particle after the statistical decay process is very 
sensitive to this difference of the internal energies of the fragments. Naive extension of the 
non-covariant QMD with the Lorentz contraction of the initial boost, called QMD/L in 
this paper, causes the unphysical instability and the spurious excitation of the fragments. 
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This leads to an unphysical result of the a multiplicity even at E\ a b= 1 GeV/u. Though 
the QMD with the Galilei boost (QMD/G) is almost free from the spurious excitation 
of the fragments, the relativistic effects in the dynamical process change the internal en¬ 
ergies of the fragments and thus clearly appear in the final a multiplicity as a factor of 
two difference compared with the QMD/G. Therefore we have to use the Lorentz covari¬ 
ant approach to analyze the fragmentation process even at the energy around f?i a b= 1 
GeV/u as long as we are concerned with the final observables of the mass distribution, 
particularly the light fragments around A = 3 ~ 4. 

In the actual analysis of the fragmentation process by the QMD and/or RQMD in 
comparison with the data, many other ingredients of the models are responsible to the 
final mass distributions. For example, the explicit momentum-dependent force has been 
already indicated to have a very important role in this energy region |ll], ^]. The Coulomb 
force is also supposed to influence the fragmentation process, particularly the expansion 
phase of the hot and/or dense matter. These are not considered in this work. Although it 
is not clear whether the relativistic effects discussed above are always visible in the final 
mass distribution or not when we include these ingredients in the models, this relativistic 
problem always exists in a fundamental point of the formulation. In order to clarify the 
effects of these other ingredients, it is also necessary to use the Lorentz covariant approach 
at least above U lab = 1 GeV/u. 
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Figure captions 


Fig. 1 : The impact-parameter dependence of the multiplicities of the a particle in 
Ca + Ca collisions at E\^= 50 MeV/u. The full squares and open circles show the results 
of the RQMD and QMD/G approaches, respectively. 


Fig. 2 : The impact-parameter dependence of the multiplicities of the proton (upper) 
and a particle (lower) in Ca + Ca collisions at F iab = 1.05 GeV/u before the statistical 
decay (left-hand side) and after that (right-hand side). The full squares, open circles, and 
diamonds show the results of the RQMD, QMD/G, and QMD/L, respectively. 


Fig. 3 : The total cross-section of fragments with the mass-number A. The full squares 
and open circles show the results of the RQMD and QMD/G approaches, respectively. 
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